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Abstract 

We analyze the 1258-day Super-Kamiokande day and night solar neutrino en- 
ergy spectra with various definitions. The best-fit lies in the LMA region 
at (Am^, tan'^^) = (5.01 x 10~^ eV^, 0.60), independently of whether system- 
atic errors are included in the x^-definition. We compare the exclusion and 
allowed regions from the different definitions and choose the most suitable 
definition to predict the regions from SNO at the end of three years of data 
accumulation. We first work under the assumption that Super-Kamiokande 
sees a flux-suppressed flat energy spectrum. Then, we consider the possibility 
of each one of the three MSW regions being the solution to the solar neutrino 
problem. We find that the exclusion and allowed regions for the flat spectrum 
hypothesis and the LMA and LOW solutions are alike. In three years, we ex- 
pect SNO to find very similar regions to that obtained by Super-Kamiokande. 
We evaluate whether the zenith angle distribution at SNO with optimum 
binning will add anything to the analysis of the day and night spectra; for 
comparison, we show the results of our analysis of the 1258-day zenith angle 
distribution from Super-Kamiokande, for which the best-fit parameters are 
(Am2,tan26') = (5.01 x 10"^ eV^ 0.56). 
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I. INTRODUCTION 



Neutrinos oscillate. Atmospheric neutrino experiments provide compelling evidence 
for this. The solar neutrino problem ||^ has been in existence for thirty years, long before 
the first indications of an atmospheric anomaly. Various solar neutrino experiments [[?|-[T^ 
detect a flux-deficit of 1/2 to 1/3 of the Standard Solar Model (SSM) prediction |13[. The 
deficit can be explained by invoking the neutrino oscillation hypothesis. Despite this, the 
solar neutrino problem is unsolved. Super-Kamiokande (SK) [0,0 has not found evi- 
dence for any of the three litmus tests for neutrino oscillations: energy spectrum distortion, 
zenith angle dependence of the flux (arising from the earth regeneration effect [|r^-|TC[) or 
seasonal variations of the flux. Moreover, there are three distinct robust solutions all of 
which have comparable significance levels, LMA, SMA and LOW [0,0. A VAC solution 
at Am2i ^ 10~^'^ eV^ is fragile in comparison to the three other solutions and its existence 



depends sensitively upon how much emphasis is placed on the SK data ^7j. We do not 
consider it further and focus on the MSW solutions. Presently, data from SK favors the 
LMA solution [jl2| . The KamLAND reactor neutrino experiment []19| will establish once and 



for all whether or not the LMA solution is correct, independent of the solar neutrino flux. 
If it is, within three years we will know Am^i and sin^2 9i2 to an accuracy of ±10% and 
±0.1, respectively |2I],2T|. The SNO experiment |^|, too, is expected to make significant in- 



roads towards the resolution of the solar neutrino puzzle p3|-p5[], especially through neutral 
current measurements. 

SNO is a Cerenkov detector with 1000 tons of heavy water as its detection medium. Its 
central objective is to test whether electron neutrinos produced in the Sun oscillate into 
active or sterile neutrinos. This can be accomplished by the simultaneous measurement of 
the rates of the charged current (CC) and neutral current (NC) reactions 

u^ + d^p + p + e- (CC) (1) 
u^ + d^p + n + u', (NC) (2) 

where denotes any of the active flavors. The NC reaction measures the total flux of 
active neutrinos which is the same as the flux produced in the Sun if there are only 
active-active neutrino oscillations. Thus, a cross section-normalized ratio NC/CC ~ 2.5 
indicates the oscillation z/g — > v^/vt- If oscillations into sterile neutrinos occur, both the CC 
and NC rates will be suppressed giving NC/CC of unity thereby signaling the existence of 
sterile neutrinos. Since the energy threshold is expected to be about 5 MeV for the CC 
reaction and 2.2 MeV for the NC reaction, only ®B and hep neutrinos will contribute to the 
SNO event rates. 

Because the NC reaction is unique to SNO, a number of studies have been devoted to 
its exploitation. In this work we undertake an analysis of what the CC rate measurements 
at SNO by themselves can and cannot tell us. They can provide valuable information for 
active-active oscillations but are not as sensitive to active-sterile oscillations. The recoil 
energy spectrum of CC events and their zenith angle distribution can in principle eliminate 
two of the three globally allowed regions in oscillation parameter space, and also measure 
the oscillation parameters. 

The stringent CHOOZ limit [^ (see also The Palo Verde Experiment |2^) of 
sin^2 0i3 < 0.1 (at the 95% confidence level), approximately decouples solar neutrino os- 
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cillations from atmospheric neutrino oscillations. For small values of ^13, provided Arn^i <^ 
Am|2 , the three-flavor survival probability P3 is related to the two-flavor survival probability 
P2 by (see the first paper of Ref. |T8|), 



P^-cos2e^^P2 0.95P2<P3<A, (3) 

where the inequality arises from the CHOOZ limit on ^13. Thus, even when the limit is 
saturated, the two-neutrino analysis represents a very good approximation to the three- 
neutrino analysis. Our analysis is performed in the two active neutrino framework. 

In Section II we briefly describe how the electron recoil energy spectra expected in the 
daytime and nighttime are calculated. We will collectively call both spectra the "D&N 
spectra". In Section III we analyze the D&N data from SK (1258 effective days) with 
different definitions and find the optimum definition for the analyses of the simulated SNO 
data. In Section IV we describe our simulation of the SNO experiment and the subsequent 
data analysis. In Section V we critically examine if the zenith angle distribution at SNO 
adds to what can be learned from the D&N spectra. We compare our expectations for SNO 
with an analysis of the zenith-angle distribution at Super-Kamiokande. We summarize our 
results in Section VI. 



II. DAY AND NIGHT RECOIL ELECTRON ENERGY SPECTRA 

The SNO CC data will provide an accurate determination of the shape of the energy spec- 
trum from neutrinos. Information about the oscillation parameters will be embedded in 
the overall suppression of the CC rate relative to that of the SSM and in the distortion of the 
shape of the energy spectrum. The low Q-value of the CC reaction (1.442 MeV) makes this 
process well-suited to obtaining a spectrum with high energy resolution because most of the 
energy of the incoming neutrino is carried away by the outgoing electron {0 < < E,^ — Q) . 
Since SNO is a real-time experiment, it is capable of studying the effect of the Earth on 
neutrinos that pass through it en route to the detector. A nadir angle, 9z, is defined as the 
angle between the negative z axis of the coordinate system at the detector and the direc- 
tion of the Sun. With this definition, cos6'^ < during the day and cos 6*^ > at night. 
Conventionally, 6z is called the zenith angle although it is actually the complement of the 
zenith angle. The relative amount of time the detector is exposed to the Sun at a particular 
zenith angle is given by the zenith- angle exposure function |28 . 

Two electron energy spectra can be measured, one each for neutrinos detected in the 
daytime and nighttime. Each electron spectrum is divided into 19 bins, every 0.5 MeV from 
the kinetic energy threshold Tth = 5 MeVQ to 14 MeV and a last bin that includes all events 
with energies from 14 MeV to 20 MeV. The expectation in each day/night bin defined by 

i?f =^ dE, (^^b{E,) + 1.8416 X 10-=^ ^hep{E,)^ P'''^(^.) cJcc{E,, AT,) (4) 



Note that the target threshold for the CC reaction is 5 MeV kinetic energy, not total energy. We 
thank E. Beier for emphasizing this point. 
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where $b and ^hep are the normahzed energy spectra of the and hep neutrinos respec- 
tively. For the undistorted spectrum shape of the neutrinos, we have adopted the result 
based on a measurement of the /5-delayed a spectrum from the decay of ®B , with this 
spectrum normalized to the flux of BPB2000. The factor 1.8416 x 10"'^ in (P is the rela- 
tive total flux of hep neutrinos to ^B neutrinos in the SSM (BPB2000). This factor is 4.5 
times larger than that of BBP98 for two reasons: (i) A recent calculation of the hep 
neutrino flux updates the BBP98 value by a factor of 4.4. (ii) In BPB2000, the ^B 
flux is 5.05 X 10^ (Itoje) cm-^ s"^ versus 5.15 x 10*^ (lto:}4) cm'^ s"^ of BBP98. The overall 
normalization A/" yields the expected number of events in the absence of oscillations if the 
Uf. survival probability at the detector, P^'^[Ey), is unity. If oscillations occur |T5|, 



(5) 
(6) 



cos2 

where Pq is the probability that a neutrino leaves the Sun as Ue, given by the well-known 



Parke formula |32| 



= l+{l-Pc^ cos2 9 cos2 9i 



(7) 



where 6^ is the mixing angle in matter at the point of neutrino production in the Sun. It is 
given by 



tan2 



tan2 6' 



1 - 



(8) 



Am*^cos2 d 



Here, is the electron density in the Sun at the creation point of the neutrino. An analytic 
expression for the crossing probability, Pc, which is a measure of the non-adiabaticity of the 



transitions is 33 



exp(-f7F)-exp(-f7-J 



1 -exp( - f 



(9) 



where 7 characterizes the adiabaticity of the resonance and F = \ — tan^^ for the exponen- 
tially varying matter density in the Sun. The adiabaticity parameter is given by [^,0 
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Am^sin^2 Q 



2E, cos2 e\N J Ne\R 



(10) 



with \Ne/Ne\R evaluated at the resonance. Finally, {P^) is the time-averaged probability 
of the transition V2 — *■ ^'e due to the effect of Earth matter^ We assume the Preliminary 



^We follow the usual conventions that vi and 1^2 are the mass eigenstates with masses mi and m2, 



the mass-squared difference, Am 



ra\, is positive and 9 is the vacuum mixing angle which 



can take values between and '7r/2, thereby accommodating an inverted mass hierarchy |42]. 
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Reference Earth Model for the Earth's electron density. The reduced cross section for 
producing an electron with measured kinetic energy in the interval ATj is 



acciK^ AT.) = dT r- dT' ^(E., T')R{T^ T') , 



(11) 



where dacc/dT' is the differential cross section for the CC reaction (from [Q) with T' being 
the actual kinetic energy of the electron. T^^j. is the kinematic limit Ey — Q. R{T, T') is the 
energy resolution function that describes the distribution of the measured energy T about 
the actual energy T' and is given by |36] 



R{T, T') 



1 



w 



V2 7rT 



:exp 



{T -T 
2^2 2-/ 



(12) 



with w = 0.348 MeV. 



III. USING THE SIMILARITIES OF SUPER-KAMIOKANDE AND SNO 

Super-Kamiokande and SNO are fairly similar experiments insofar as CC measurements 
are concerned. Both are high-statistics real-time electronic experiments using Cerenkov light 
detection and both are sensitive to only and hep neutrinos because their energy thresholds 
are almost the same. SK uses H2O as its detection medium while SNO uses D2O. Thus, 
SK detects solar neutrinos via the elastic scattering (ES) reaction 

iy, + e' + e- (ES) , (13) 

which only makes a minor contribution to the rate at SNO. However, since the neutrino 
source and the principle of neutrino detection are the same in both experiments, it is reason- 
able to expect the two experiments to yield equivalent flux measurements. This equivalence 
has been exploited to devise ways to predict the NC rate at SK once SNO has CC rate 
results ||3^ and to predict the energy spectrum at SNO from that measured by SK |3^ 



We assume that the equivalence of SK and SNO is sufficiently robust that the best 
definition for SK will also be the best for SNO. The calculation for the ES rate at SK is 
similar to the CC rate described above for SNO except for the following alterations: 

(i) P occ is replced by Pde + (1 — -P) ct/x, where and are the z/g — e and i/^ — e ES 
cross sections [^, respectively. 



(ii) The bins are defined in terms of the total electron energy, since SK reports its data in 
terms of the reconstructed total energy of the recoil electron, with threshold Ey = h 
MeV. 

(iii) In Eq. (0), = 0.47 MeV |n|. 



SK has reported results from 1258 days of data-taking [|ri|,|T^] as ratios with respect 
to the first version of BPB2000 in which the ^B fiux is 5.15 x 10^ (ll[!:?^) cm'^ s'^ |9 . 
We call this SSM'. A recently revised version of BPB2000 gives the '^B fiux as 
5.05 X 10^ (l^aie) cm~2s^^ ||13[. Consequently, we modify the SK data accordingly, by mul- 
tiplying the central value and statistical error in each bin by the ratio 5.15/5.05. The 
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systematic errors are conveniently given as percentages and do not need modification. The 
measured flux suppression is [lTTlp!2 



which relative to the SSM is 



D^^^SK _ „ +0.017 

SSM' - U-4^^-0.015 , 



P^^^SK _ „ . c-Q+0.018 

SSM - ^-^^^-0.016 • 



(14) 



(15) 



SK has presented results using several different definitions. In their latest flux- 



independent analysis of the D&N spectra, they used [[12 



^^^(Am^tan^ 



38 

E 



' ^uncorr'j 2 



(16) 



where the flux measured by SK in the i*^ bin is 0™^^"*, the expected flux without oscillations is 
(pSSM ^YiQ expected flux with oscillations is 0°**^ = 0°*'^(Am^, tan^^). The uncertainties 
^stat^ ^uncorr ^corr statistical, uucorrelatcd and correlated uncertainties in the 

^th Y)^^^ respectively. The correlated errors include the experimental uncertainties in the 



determination of the spectrum and the theoretical uncertainties in the calculation 
of th expected energy spectrum The functions, /i(/5), parameterize the correlated 



uncertainty in the shape of the spectra; /? is a free shift factor of the correlated error and 
a is a free parameter that normalizes the measured flux relative to the expected flux. The 
sum runs over 38 energy bins (19 day bins + 19 night bins). 

We have taken the SK 95% C. L. exclusion region from Ref. ||12| for our later comparison 
with regions obtainable with alternative deflnitions. It is the hatched region enclosed by 
the dotted line in Fig. Note that this exclusion region corresponds to data relative to 
SSM'. The dark shaded areas are the allowed regions at 99% C. L. from a global analysis 
with free ®B and hep fluxes and the SSM (not SSM') [|l^]. The LOW solution is allowed only 
at the 99% C. L.. The analysis includes the D&N spectra from the 1117-day event sample 
but not the total rate, since that information is contained in the energy spectra. Ref. [l^ 
used the shape of the undistorted spectrum of *B neutrinos from Ref. |^ . We emphasize 
that although the SMA region found from a combined flt of flux measurements is excluded 
at the 95% C. L. ]g, the SMA region from the global flt of Ref. |3 is not. 

Another suitable deflnition of x^, similar to one used by SK in earlier analyses (of 825 
effective days of data), is 



X^(Am^, tan^ 



38 

E 



^oscl^SSM^ 



(17) 



Again, a is a free flux normalization factor and (3 constrains the variation of correlated 
systematic errors. Performing a analysis of the 1258 day data with this deflnition gives 
the 95% C. L. exclusion (hatched) region (x^ > 52.19 for 37 d.o.f.), outlined by the dashed 
line in Fig. ^ The hep contribution to the neutrino flux is left unconstrained in the SK 
analysis while we have flxed the ratio between the hep and ^B fluxes. Thus, for a test of the 
MSW hypothesis, the SK analysis has 36 degrees of freedom while our analysis has 37. Also, 
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keeping in mind that the regions enclosed by the dotted and dashed hnes correspond to two 
different reference solar models and definitions, it is noteworthy that the general shapes 
of the regions closely resemble each other, although they differ in size as a consequence of 
the term J2 {(^/(^T^^Y x\k- Even if j3 is very small, (3 ~ 0.001, the contribution from this 
term can significantly increase the value of thereby permitting a larger exclusion region^. 
On the other hand, by including the correlated errors in as in (|T7|) , their effect is greatly 
diminished. It is evident that the spectral distortion functions fi{l3) play an important role 
in defining an efficient function. To include the possibility of negative Am?, we have also 
plotted the same region (dashed line) with tan^6' as the abscissa 

In the approximation that systematic errors can be neglected, both the above defini- 
tions lead to 



38 r 

xLi(Am2,tan2^) = ^ 

i=l 



j^meas j ^SSM 



a(f)[ 



osc/^SSM 



stat 



(18) 



The resulting 95% C. L. exclusion region (hatched and enclosed by the solid line) is shown in 
Fig. ^ Again, we note the remarkable similarity of the shapes of the exclusion regions from 
the three analyses. Dropping systematic errors leads to a region more similar in size to that 
obtained by the SK collaboration than that obtained by using Eq. (|T7p. On this basis, we 
hereafter assume that we can safely ignore all systematic errors when making projections for 
SNO with simulated data. If SK is a reasonable guide, we will err on the conservative side. 
It may be counterintuitive that the exclusion region using x^sk larger than that using x^stat 
because one expects more errors to lead to less confidence and therefore a smaller exclusion 
region. However, as explained earlier, the correlated errors in Eq. |1^ are responsible for this. 
The region using Eq. (|l^) is smaller than that of x^stati agreement with expectations. 

So far we have only considered flux-independent exclusion plots. If instead, flux- 
dependent allowed regions are sought, one needs to add another term to the x^ definitions 
considered, 



X^(Am^, tan^ 



^^(Am^tan^^ 



+ 



a 



On 



(19) 



where =l'lo;i6SSM (or SSM') is the theoretical uncertainty in the *B flux. In our analysis 
we symmetrize this value to cTq = ±0.18 SSM. The 95% C. L. allowed regions (Ax^ < 5.99 
for two oscillation parameters), are superimposed on the exclusion plots in Fig. The 
allowed region from SK's analysis is not shown in Ref. [|l^]; we have taken it from Ref. |^ 



With the 
.tan^^) = 
18|) for 36 



It is evident that the allowed regions are alike with minor differences in size. 
X^-definitions of Eqs. (0) and (pISD, we find the same best-fit parameters, (Am^ 
(5.01 X 10"^ eV^ 0.60) with x^ = 30.6 using Eq. (P) and x^ = 32.4 using Eq. 
degrees of freedom. 

In Fig. ^ the crosshairs represent the best-fit parameters which are very close to those 
presented by SK from an analysis of the data from 1117 days in which they include the flux 



^If we set fi{(3) = 1, thereby making the flux normalization the same in all bins, and /? 
we improve the agreement with the SK region significantly. 



0.001, 
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constraint SK does not report the best-fit point from a flux-dependent analysis of the 
1258-day D&N spectra. From a flux-independent analysis they find that the minimum ^ 
value lies in the VAC region [Q, which we have not considered in our analysis. However, 
from a flux-dependent analysis of their zenith spectrum, they find the minimum x 
VAC region, with some points in the LMA region with similar values |1 
(Am^,tan^^) = (7 x 10~^ eV^, 0.47) is one such point, which is close to our best-fit param- 
eters. A cautionary note when interpreting Fig. |I] is that the dark-shaded flux-independent 
globally allowed regions of Ref. [1^ are not directly comparable to our flux- dependent al- 



^ in the 
For example, 



lowed regions. Our motivation for superimposing the flux- independent allowed regions is to 
facilitate a comparison with the flux-independent exclusion regions. 



IV. DATA SIMULATION AND ANALYSIS 

If the SSM flux is correct, then in the absence of oscillations SNO should detect about 
9250 events per year. If instead the SSM flux is wrong and no oscillations occur, the flux 
suppression expected at SNO is the same as that seen by SK, Eq. (|T3p. We first simulate 
data assuming this pessimistic scenario and predict the exclusion and allowed regions that 
we can expect SNO to present three years from now. Then we turn to the more reasonable 
explanation in which neutrino oscillations do occur. 

We simulate data for the best-fit oscillation parameters in each of the three allowed MSW 



solutions (from the global analysis of Ref. |T^) and display the corresponding exclusion and 



allowed regions. Exclusion regions are meaningful only when the data points are normally 
distributed and that for at least some region of the parameter space of interest, x^/d.o.f ~ 1. 
We enforce this by simulating data for which this is true at the input oscillation parameters. 
The normalization constant M in Eq. (^) is set by the stipulation that we are considering 
three years of data accumulation. 

Figure | shows the expected spectra (as a ratio with respect to the SSM) for typical LMA 
(solid histogram), SMA (dashed histogram) and LOW (dotted histogram) solutions. The 
data points here are simulated for the LMA solution. For direct comparison of the spectral 
shapes, these spectra are normalized to the flux corresponding to the simulated data. 

For the sake of specificity, we define as 



X^(Am^tan2^)=$: 



38 r J^simulated j ^SSM _ ^ /^o«c j ^SSM -| 2 



1=1 



^stat 



(20) 



where af""^ = y'/?|™"'»*e(i^^5S'M ^^^^ number of simulated events, /?|«™"'«*ec(^ 
ATj is obtained by randomly choosing a point from a Gaussian distribution centered at 
the theoretical value and of width equal to the square root of the theoretical value. This 
definition is the same as that of Eq. (|19D except that it is expressed in terms of the number 
of events rather than the flux. 

To be conservative, we only show 99% C. L. exclusion (x^ > 59.89 for 37 d.o.f.), and 
allowed regions (Ax^ < 9.21), resulting from the simulated SNO data. Figure || shows the 
expected regions for the flat spectrum hypothesis (with the same flux suppression as seen 
by SK) and three MSW solutions. The stars and crosshairs mark the theoretical inputs and 
best-fit points, respectively. The SMA expectation is characteristic and easily identifiable. 
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Dataset 



Day 


-1 


< 




< 





Nl 





< 


cosOz 


< 


0.173 


N2 


0.173 


< 


cosOz 


< 


0.5 


N3 


0.5 


< 


cos9z 


< 


0.707 


N4 


0.707 


< 


cos9z 


< 


0.83 


N5 


0.83 


< 


cos 9z 


< 


0.92 



TABLE I. The definitions of the night bins in terms of the nadir of the Sun, 9z- Note that the 
"core bin", N5, does not contain zeniths beyond 0.92 because the latitude of the detector restricts 
its range. 



For all other possibilities, the plots bear a striking semblance to each other and to the SK 
results. The LMA and LOW solutions persist simultaneously for all but the SMA solution. 
The SMA region is excluded to a large extent and with a more efficient definition, the 
entire SMA region could be excluded. 



V. ZENITH ANGLE DISTRIBUTION 

In principle the information contained in the zenith angle distribution is contained in the 
D&N spectra. The SK exclusion regions obtained by them in separate analyses of the D&N 
spectra and the zenith angle distribution (with energy subdivisions) bear testimony to this 
expectation; the differences are small [0. This is despite the fact that each zenith angle 
bin is split into several energy bins, thereby potentially maximizing the resolution available 
to SK. It has been advocated that an appropriate choice of binning might make it possible 
for SNO to not only identify which solution is correct but also to determine the oscillation 
parameters . By performing a complete analysis, we now assess the extent to which this 



claim can be validated. 

The choice of night bins of equal size in cos 6z is democratic, but does not take advantage 
of any distinctive features of the distributions of the different solutions. In the context of 
SNO, the qualitative behavior of the distributions for the LMA, SMA and LOW solutions 
was studied in detail in Ref. 0]. It was found that with a suitable choice of binning, 



any smearing of peculiarities intrinsic to the region of parameter space can be avoided. 
The events were binned in a manner that leads to a characteristic distribution for the LOW 
solution because in the SMA and LMA regions, the cos ^^-dependence is rather weak leading 
to a more or less flat distribution. This remark is pertinent, since as we have seen, it is 
difficult to differentiate the LMA and LOW solutions from each other. If the SMA solution 
is the correct one, the strong spectral distortion will easily make it stand apart. As deflned 
in Table I, there is one day bin and flve non-uniform night bins. The "core bin", N5, at 
SNO {{cos 9z)max = 0.92) is smaller than that at SK {{cos 9z)max = 0.975) because SNO's 
higher latitude restricts the 9z range. Thus, a smaller number of solar neutrinos that pass 
through the Earth's core are incident at SNO than at SK. 

Since we know the zenith-angle exposure function ||2^, the SSM prediction for the number 



of events in each zenith angle bin, RfssMi ^-nd the prediction with oscillations, 
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be calculated. For a given set of oscillation parameters, we want to generate zenith angle 
distributions that have the same number of events as the simulated energy spectra of the 
previous section. The number of events in the day bin is simply the sum of all the events 
in the day spectrum. For the night bins, we use Rf osc central value of a Gaussian 

distribution and simulate the number of events in each night bin. Note that the number 
of events at night represented by this distribution does not coincide with that of the night 
spectrum. The nighttime distribution is renormalized to yield the number of events in the 
night spectrum. Now the simulated energy spectra and zenith angle distribution reflect the 
same data. To match the number of simulated events in the zenith angle distribution and 
the D&N spectra, two normalizations are needed, one each for the daytime and nighttime 
events. The shape of the theoretical expectation, however involves only one normalization, 
namely, the total number of simulated events. 

We perform fiux-independent analyses (for the same datasets used to find the regions of 
Fig. ^ with the simple function. 



X^(Am^, tan 



2 i.^^2c^Y ^ ^ 

i=l 



Rf, simulated/ Rf, SSM ^ Rf,osc/Rf,SSM 



cyf 



(21) 



where of = y Rf simulated/ Rf.ssM- From Fig^, it is evident that the zenith angle distributions 
of the various solutions lead to similar 99% C. L. exclusion regions (x^ > 15.09 for 5 d.o.f.). 
We underscore the fact that the regions corresponding to the flat spectrum, LMA and 
LOW datasets exclude part of the SMA region, and that the region found with the SMA 
dataset excludes a large part of the LMA and LOW solutions. This is consistent with 
Fig. ^. Additionally, the 99% C. L allowed regions have the same shapes for the LMA and 
LOW solutions. The flux constraint is included to find the allowed regions. The stars and 
crosshairs mark the theoretical inputs and best-fit points, respectively. 

Figure |^ shows the 99% C. L. exclusion (x^ > 16.81 for 6 d.o.f.), and allowed regions 
from the zenith angle distribution from 1258 days of data at SK. In making these regions, 
we have employed SK's binning which is different from our choice for SNO. The statistical 
and systematic errors are added in quadrature. The crosshairs mark the best-fit parameters, 
(Am^tan^^) = (5.01 x 10"^ eV^ 0.56) with = 5.15. The exclusion regions of Figs. ^|-|^ 
are similar to the 99% C. L. exclusion region reported by SK (with a 504-day dataset) 
with a similar definition No parts of the globally allowed regions are excluded. Of 
course, subdividing each zenith angle bin into energy bins (as done by SK in their latest 
analysis ||12|) will greatly improve the sensitivity of the analysis, but this is equivalent to 
using D&N energy spectra. 

The day-night variation embodied in the zenith angle distribution can also be presented 
in terms of a day-night asymmetry defined by 

N - D pN _ pD 
= 2 = 2 — , (22) 

where D and N are the total number of events detected during the days and nights, re- 
spectively. The approximate ranges of Adn are (0.005,0.1) in the LMA and LOW regions 
and (—0.01, 0.05) in the SMA region. Note that in parts of the SMA region, A^n < 0, thus 
uniquely identifying the SMA solution. However, Aon > —0.01 and an identification of such 
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a small deviation from zero will be difficult. Semi- analytic approximations for Adn have 



been derived for the three allowed regions [^. These expressions can be used to provide 



insight into the orientations of the zenith-angle allowed regions. Iso-asymmetry lines in A^itv 
that pass through all three regions must have A^jy > and are given by 

Am^ (eV^) ^ 3 x IQ-^ (LMA) , (23) 

AoN^sm^^e (SMA), (24) 

Am^ (eV^) ^ 2.5 x 10"^ (LOW) . (25) 

sin 2 9 



The relations in Eqs. (P3H25|) have a wider domain of applicability than indicated. For 
example, Eq. (|25| ) is applicable in the range 10~^ < Am^ ^ 3 x 10~^. This explains why 
the 99% C. L. allowed regions almost connect the SMA and LOW solutions. 



VI. SUMMARY 

We have analyzed the 1258-day day and night energy spectra presented by Super- 
Kamiokande using definitions that account for systematic errors in different ways. The 
best-fit lies in the LMA region at (Am^,tan^6') = (5.01 x 10~^ eV^, 0.60), independently 
of whether systematic errors are included in the x^-definition. We have shown that these 
approaches lead to exclusion and allowed regions of different sizes, but the general areas of 
the regions remain unchanged even if systematic errors are neglected (see Fig. [^). Using 
Super-Kamiokande as our reference, we then draw conclusions for SNO based on analyses 
that incorporates only statistical errors. We assume the optimistic electron kinetic energy 
threshold of 5 MeV and three years of accumulated data. 

If the SMA solution is correct, the day and night spectra will show sufficiently strong 
distortions to distinguish it from the other solutions (Fig. 0(d)). For a flux-suppressed fiat 
spectrum or the LMA and LOW solutions, the regions are similar enough to not provide any 
constraint beyond the exclusion of most of the the SMA solution (Fig. ^(a-c)); the LMA 
and LOW solutions are indistinguishable. However, KamLAND will certainly help in 
this regard by either ruling out the LMA solution or pinning down the LMA oscillation 
parameters ||20|. The zenith angle distribution in itself will not add anything to what can 
be obtained from the day and night spectra unless each zenith angle bin is subdivided into 
energy bins. 

The 99% C. L. exclusion regions in Fig. ^bear a striking similarity to that obtained from 
the 1258-day zenith angle distribution at Super-Kamiokande (Fig. §), for which the best-fit 
parameters are (Am^,tan^^) = (5.01 x 10~^ eV^, 0.56). We expect results from charged- 
current measurements at SNO to be similar to that of Super-Kamiokande, thus providing 
an important check of the Super-Kamiokande conclusions. Needless to say, neutral current 
data from SNO will provide a crucial test of the existence of sterile neutrinos, and if solar 
neutrinos do not oscillate to sterile neutrinos, SNO will measure the flux produced in the 
Sun. 
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SK D&N Spectra {u^ — - ^^.t) 



1-3 





tan^e 

FIG. 1. The exclusion (hatched) and allowed (lightly shaded) regions at 95% C. L. obtained 
from the 1258 day SK D&N energy spectra using three different definitions. The regions enclosed 
by the dotted (taken from Ref. |l^), dashed and solid lines result from the use of Eqs. (p^),(p^ 
and ([I^), respectively. Equation ( p!^ ) has no contribution from systematic errors. The flux 
constraint is imposed to find the allowed regions. The crosshairs mark the best-fit parameters, 
(Am^,tan^6') = (5.01 x 10~^ eV^, 0.60) using Eqs. (0) and (p^). The dark shaded regions are the 
global solutions (with 1117 SK days) at 99% C. L. with free and hep fiuxes found in Ref. |17|. 
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FIG. 2. The expected e~ D&N energy spectra at SNO (as a ratio with respect to the SSM) 
with three years of accumulated data for typical LMA (solid histogram) , SMA (dashed histogram) 
and LOW (dotted histogram) solutions. The data points are simulated for the LMA solution. The 
spectra are normalized to the flux reflected by the data. The last bin includes all energies from 
14 MeV to 20 MeV. 
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FIG. 3. Expected 99% C. L. exclusion (hatched) and allowed (lightly shaded) regions from the 
D&N spectra at SNO for (a) a flux-suppressed flat spectrum (b) LMA solution, (c) LOW solution 
and (d) SMA solution with three years of accumulated data. The dark shaded regions are the 
global solutions at 99% C. L. with free '^B and hep fluxes. The stars and crosshairs mark the 
theoretical inputs and best-flt points, respectively. 
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FIG. 4. Expected 99% C. L. exclusion (hatched) and allowed (lightly shaded) regions from 
zenith angle distributions at SNO corresponding to the simulated datasets of Fig. ^. The stars and 
crosshairs mark the theoretical inputs and best-fit points, respectively. 
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SK Zenith. Angle Distribution (Ug — -^/i.t) 




tan^e 

FIG. 5. The exclusion (hatched) and ahowed (hghtly shaded) regions at 99% C. L. from 
the 1258-day SK zenith angle distribution. The crosshairs mark the best-fit parameters, 
(Am2,tan26l) = (5.01 x 10"^ eV^ 0.56). 
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